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ABSTRACT 


A  multiple-wavelength  light  transmittance  measurement 
system  previously  used  in  a  laboratory  environment  to  study 
particles  in  solid  rocket  propellant  exhaust  plumes  was 
modified  for  use  in  the  field,  where  high  levels  of  vibration 
can  degrade  the  accuracy  of  data.  The  system  was  converted 
from  a  single  light  beam  configuration  to  a  dual  beam 
configuration  which  was  capable  of  obtaining  a  complete  set  of 
1024  reference  and  scene  measurements  in  10  ms.  Modifications 
included  designing,  building  and  testing  a  new  analog-to- 
digital  data  converter  trigger  circuit,  and  a  rotating-wheel 
light  chopper.  Optical  components  including  beam  splitters, 
lenses  and  a  fiber  optic  cable  were  installed,  and  existing 
data  collection  system  software  was  modified.  The  new  system 
was  tested  by  measuring  soot  from  an  oxy-acetylene  torch  to 
prove  the  design  concept.  Test  results  and  system  performance 
were  documented.  Recommendations  for  further -modifications, 
improvements  and  applications  are  presented. 
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I .  INTRODUCTION 


During  a  recent  war,  a  need  for  exhaust  plume  signature 
prediction  and  identification  for  liquid  propellant  rocket 
engines  was  realized.  Since  radiation  emissions  are  greatly 
effected  by  particles  in  the  plume,  knowledge  of  their  size, 
optical  properties,  and  concentration  provides  valuable 
information  which  is  needed  to  predict  the  exhaust  signature. 
Liquid  rocket  propellants,  unlike  solid  rocket  propellants,  do 
not  normally  contain  any  added  metals  which  can  enhance  the 
signature  of  the  plume.  However,  the  exhaust  can  contain 
significant  amounts  of  carbon  particles,  which  result  from  the 
combustion  of  hydrocarbons.  Soot  can  produce  increased  visible 
and  infrared  signatures. 

Several  methods  for  measuring  particles  in  exhaust  plumes 
exist.  The  most  common  procedures  are  particle  sample 
withdrawal  followed  by  optical  measurements  or  in  situ  methods 
based  on  the  transmission  and/or  scattering  of  iight  [Ref.  1]. 
Soot  particles  produced  by  liquid  propellants  are 
predominately  in  the  submicron  range  [Ref.  2,3].  Measurement 
of  these  small  particles  can  be  accomplished  using  multiple- 
wavelength  light  transmission  techniques  [Ref  4]. 

Recent  research  at  the  Naval  Postgraduate  School  by  Vaughn 
[Ref.  1]  and  Kim  [Ref.  5]  successfully  used  multiple- 
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wavelength  light  transmission  measurements  for  submicron 
aluminum  oxide  particles  in  the  plumes  of  solid  propellant 
rockets.  This  technique  measured  the  transmittance  at  six 
wavelengths .  The  measured  ln-transmittance  ratios  were 
compared  to  the  extinction  coefficient  ratios  obtained  from  a 
Mie-based  computer  code.  The  Beer-Lambert  law  for  the 
exponential  decay  of  light  intensity  through  a  uniform 
particle  field  could  then  be  used  to  obtain  the  particle 
concentration  [Ref.  4]. 

The  apparatus  used  in  these  investigations,  while  accurate 
for  small  rocket  motors  in  a  controlled  laboratory 
environment,  would  probably  be  greatly  effected  by  vibrations 
caused  by  larger  motors  [Ref.  2].  The  apparatus,  as 
originally  configured,  collected  all  reference  measurements 
for  transmittance  calculations  prior  to  rocket  motor  ignition, 
which  were  then  averaged  to  get  a  mean  value .  Light 
measurements  through  the  plume  were  then  collected  after  motor 
ignition,  which  were  also  averaged.  Transmittance  was  obtained 
from  the  averaged  values  [Ref.  1,5]. 

Vibrations  from  a  large  motor  firing  can  cause  variations 
in  the  intensity  of  the  light  source  and/or  the  alignment  of 
the  optical  components,  thus  invalidating  the  averaged  pre-run 
reference  measurements  [Ref.  2].  One  way  to  overcome  this 
limitation  is  to  collect  reference  measurements  simultaneously 
with  plume  measurements.  Transmittance  can  then  be  calculated 
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at  any  instance  of  time,  thus  negating  the  effects  of 
vibration  and  other  variations  during  the  burn  time.  These 
individual  transmittances  could  then  be  averaged,  if  desired. 

Another  limitation  of  the  light  transmission  and 
collection  equipment  was  that  it  had  only  been  tested  across 
short  distances  (less  than  two  feet)  [Ref.  1,5 J.  Use  of  the 
apparatus  on  full  scale  rocket  motors  would  require 
transmission  lengths  of  fifteen  feet  or  more. 

The  purposes  of  this  investigation  were  to  (1)  modify  the 
multiple-wavelength  light  transmission  apparatus  and 
associated  data  collection  equipment  used  by  Vaughn  [Ref.  1] 
to  obtain  near-simultaneous  reference  and  plume  measurements, 
(2)  extend  the  light  transmission  distance  and  (3)  validate 
the  modifications  by  measuring  soot  particles  from  a  test 
source. 
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II.  THEORETICAL  BACKGROUND 


A.  INTRODUCTION 

Multiple  wavelength  light  transmission  measurement 
techniques  have  been  used  extensively  in  recent  years  to  study 
the  nature  of  small  particles.  The  theory  behind  these 
techniques  is  discussed  in  detail  by  Few  [Ref.  2,6], 
Cashdollar  [Ref.  4]  and  Vaughn  [Ref.  1]. 

Recent  particle  sizing  techniques  using  light  transmission 
measurements  are  varied  applications  of  the  Mie  theory  for  the 
light  scattered  by  particles  [Ref.  1,2].  For  particles  much 
smaller  than  the  wavelengths  of  the  incident  light  the 
Rayleigh  approximation  is  often  used.  The  techniques  have 
been  refined  over  the  years,  but  the  basic  approach  is  that  a 
collimated  light  beam  is  transmitted  through  a  particle  field 
[Ref.  2].  Particle  size  information  is  then  provided  by  the 
wavelength  dependence  of  the  angular  scattering  and  extinction 
of  the  light  beam  [Ref.  2]. 

The  technique  used  in  this  experiment  focused  on  using 
the  extinction  of  light  of  multiple  wavelengths  to  determine 
mean  particle  size  and  the  particle  index  of  refraction. 


4 


B.  LIGHT  TRAMSMISSIOM/EZTIMCTIOH  THEORY 

The  transmission  of  light  through  a  uniform  particle  cloud 
is  given  by  the  Beer-Lambert  law  [Ref.  4]: 


T=exp ( -QAnL) =exp ( - 


2pd 


(1) 


where  T  =  fraction  of  light  transmitted, 

Q  =  dimensionless  extinction  coefficient, 

A  =  cross  sectional  area  of  a  particle, 
n  =  number  concentration  of  particles, 

L  =  path  length, 

C.  =  mass  concentration  of  particles, 
p  *  density  of  an  individual  particle, 
d  =  particle  diameter 

Mie  scattering  theory  considers  the  interaction  of  a  plane 
electromagneteic  wave  with  a  spherical  object.  Computation  of 
Q  as  a  function  of  particle  size,  light  wavelength,  and 
complex  index  of  refraction  is  possible,  and  further 
extensions  of  the  theory  have  been  applied  to  nonsperical 
particles  [Ref.  1,2,4]. 

This  extension  holds  if  the  nonspherical  particles  are 
approximated  by  spherical  particles  whose  diameters  provide 
the  same  volumes  as  the  nonspherical  particles.  [Ref.  1,2,6] 
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For  a  polydisperse  system  of  particles,  Dobbins  [discussed 


Ref.  4]  revised  the  transmission  law  to  obtain: 


,  3<5C_Lt 

r=exp(-  -  — ;  - ) 


2p  d- 


(2) 


32 


where  Q  is  the  average  extinction  coefficient  and  d32  is  the 
volume-to-surface  mean  (Sauter  mean)  particle  diameter  [Ref. 
1/4]. 

Taking  the  natural  logarithm  of  equation  ( 2 )  yields : 


inr'5(-^L) 

2pa32 


(3) 


A  multiple-wavelength  light  beam  has  identical  L,  p,  C,,  and 
d32  for  each  wavelength  through  the  plume.  Therefore,  the 
ratio  of  the  ln-transmissions  at  any  two  wavelengths  is  equal 
to  the  ratio  of  the  calculated  mean  extinction  coefficients 
for  the  same  wavelengths  [Ref.  1,4]: 

lnTUj) 

lnT(X.)  _  d,,) 

J  experimental  32'  theoretical 
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Cashdollar's  computer  code  [discussed  in  Ref.  1,4],  utilizing 
Mie  theory,  was  used  to  generate  the  values  of 

Transmittances  at  the  various  wavelengths  are  determined 
experimentally.  Previous  transmittance  measurements  by  Vaughn 
[Ref.  1]  and  Kim  [Ref.  5]  used  a  single  light  beam,  as 
discussed  in  Chapter  I.  In  a  dual  beam  system,  as  used  by  Few 
[discussed  in  Ref.  2],  a  calibration  measurement  of  the  ratio 
of  scene  beam  (S)  intensity  to  reference  beam  intensity  (R)  is 
conducted  prior  to  the  introduction  of  soot  particles, 
yielding  [Ref.  2]: 


Real 


ISU) 

J,U) 


(5) 


The  same  intensity  ratio  is  then  recorded  with  soot  particles 
present  in  the  scene  beam.  The  transmittance  of  light  at  each 
wavelength  is  obtained  as  the  ratio  of  the  sample  ratio  to  the 
calibration  ratio  [Ref.  2]: 


rU)  = 


( x ) 


tt) 


(6) 
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The  procedure  for  determining  soot  particle  diameter 
requires  calculation  of  Q  for  various  assumed  values  of  k,  d32, 
m=a-ib  (complex  index  of  refraction),  and  an  assumed  particle 
size  distribution.  If  the  assumed  values  are  correct, 
calculated  Q  ratios  will  match  the  measured  ln-transmission 
ratios  for  each  of  the  wavelengths  utilized.  [Ref.  1,4] 

Studies  by  Cashdollar  [Ref.  4],  Few  [Ref.  2],  and  others 
[discussed  in  Ref.  2,4]  report  that  the  complex  index  of 
refraction  of  soot  has  a  real  part  (a)  value  of  1.8  to  1.9. 
The  complex  part  (b)  has  values  ranging  from  0.3  to  0.66. 
Additionally,  the  studies  report  that  small  particles  in 
plumes  often  have  a  monomodal  log-normal  distribution.  In 
this  case  the  distribution  is  characterized  by  d32  and  o 
(standard  deviation).  [Ref.  1,4] 

Since  there  are  four  unknown  variables  (d32,  o,  a,  and  b) 
a  minimum  of  four  independent  ln-ratios  are  needed.  Increased 
numbers  of  ln-transmission  ratios  increase  the  ability  to  find 
the  best  fit  between  experiment  and  theory.  The  values  of  d32, 
o,  a,  and  b  are  found  by  finding  the  best  fit- (in  the  least- 
squares  sense)  between  the  calculated  and  measured  results. 
[Ref.  1,2] 
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III.  EXPERIMENTAL  APPARATUS 


A.  GENERAL  ARRANGEMENT 

Figure  1  shows  the  general  arrangement  of  the  light 
transmission  and  collection  equipment.  A  variable  beam 
splitter,  directly  in  front  of  the  light  source,  was  used  to 
split  the  light  into  scene  (through-plume)  and  reference  beams 
for  transmittance  measurements .  The  reference  beam  was 
redirected  through  a  chopper  and  focused  into  a  fiber  optic 
cable,  which  was  routed  around  the  particle  cloud.  The  scene 
beam  passed  directly  through  the  chopper,  a  collimating  lens, 
and  then  through  the  plume.  The  chopper  provided  alternating 
light  measurements  for  the  reference  and  scene  beams .  A 
second  beam  splitter  was  used  to  direct  the  reference  beam 
from  the  fiber  optic  cable  into  the  spectrograph. 
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Figure  1 

General  Arrangement 
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A  200  Watt  Mercury-Xenon  bulb  (Oriel  Corp.  #6291)  was  used 
as  the  light  source  in  this  experiment.  It  produced  intensity 
peaks  at  wavelengths  of  253.4,  302.2,  312.6,  334.2,  366.3, 
404.7,  435.8,  546.1,  and  577.0  nm.  The  bulb  was  mounted  in  a 
Oriel  lamp  housing  with  a  F/l  condenser  and  a  UV  grade  fused 
silica  1.5  inch  lens.  The  lamp  was  powered  by  an  Oriel  power 
supply.  [Ref.  7] 

An  Oriel  model  77400  Multispec  spectrograph  was  used  to 
collect  the  reference  beam  and  scene  beam  light  intensities. 
A  25pm  wide  input  slit  restricted  the  the  amount  of  stray 
light  entering  the  spectograph.  The  transmitted  light  was 
then  projected  onto  a  1024  element  diode  array.  Figure  2 
shows  the  optical  pattern  within  the  spectrograph. 

The  diode  array  was  an  EG&G  Reticon  G  series  solid  state 
line  scanner  with  1024  elements.  Spectral  response  of  the 
array  3  in  the  UV/IR  range  (approx.  200  -  1000  nm 
wavelength).  The  diode  array  was  controlled  by  a  Reticon 
RC100  motherboard  which  provided  all  clock,  start,  video 
amplifier,  and  blanking  requirements.  [Ref.  81 
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EXIT  slit,  or  flange 
for  oiode  array 


Figure  2 

Multispec  optical  arrangement 
{Ref.  7] 


C.  LIGHT  CHOPPER 


A  detailed  drawing  of  the  light  chopper  disk  is  provided 
in  Appendix  A.  The  two  sets  of  large  apertures  were  designed 
to  provide  alternating  reference  and  through-plume  (scene) 
light  measurements.  The  aperture  arc  length  provided  a 
transmission  time  of  slightly  more  than  the  sweep  time  of  the 
diode  array  (2.04  ms)  when  the  disk  was  rotated  at  1000  rpm. 
The  small  apertures  on  the  outer  radius  provided  a  "take  data" 
command  signal  to  the  Data  Acquisition/Control  System  and  were 
aligned  with  the  leading  edges  of  the  large  apertures.  The 
light  chopper  provided  1024  transmittance  values  each  10.0  ms. 

The  chopper  disk  was  driven  by  a  Vexta  model  FBL-215A  24 
volt  DC  speed  control  motor.  The  motor  speed  was  controlled 
by  a  precision  potentiometer.  Speed  was  monitored  by  checking 
the  frequency  of  the  "take  data"  command  from  the  chopper 
disk. 

D.  FIBER  OPTICS 

The  reference  beam  was  routed  around  the  particle  cloud 
using  an  Oriel,  600  pm  diameter,  UV/VIS  single  fiber  optic 
cable.  Cable  length  was  5  meters.  Light  was  collected  into 
the  cable  using  an  Oriel  model  77780  Fiber  Optic  Input 
Accessory  with  a  fused  silica  F/2  lens.  Output  from  the  cable 
was  directed  onto  the  second  beam  splitter  using  an  Oriel 
model  77644  Collimating  Beam  Probe  with  fused  silica  lens. 
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E.  OTHER  OPTICS 

Two  types  of  beam  splitters  were  used.  The  first  type 
was  a  rotating  disc,  variable  beam  splitter  which  allowed  for 
proportioning  appropriate  intensities  of  light  into  the  scene 
and  reference  beams.  A  fused  silica  50/50  polka-dot  beam 
splitter  was  placed  in  front  of  the  spectrogaph.  It  was  used 
to  redirect  the  reference  beam  into  che  spectrograph  entrance 
slit. 
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IV.  DATA  ACQUISITION  SYSTEM 


A.  OVERVIEW 

Data  and  clock/timing  signals  from  the  diode  array  were 
sent  to  a  Hewlett  Packard  6942A  Multiprogrammer  where  the  data 
were  converted  from  analog  to  digital  and  stored  on  memory 
cards  for  processing.  The  multiprogrammer  was  programmed  by 
a  HP  9836S  desktop  computer.  Figure  3  is  a  block  diagram 
showing  the  general  arrangement  of  the  data  acquisition 
system. 
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Figure  3 

Data  Aquisition  System 
[Adapted  from  Ref.  9:p.  40] 
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B.  MULTI PROGRAMMER 


The  HP6942A  Multiprogrammer  is  a  self-contained  system 
which  has  a  wide  variety  of  data  acquisition  applications. 
The  system  can  accomodate  up  to  sixteen  I/O  cards  which  are 
programmable  by  a  desktop  computer  through  the  HP  Interface 
bus.  The  cards  can  be  addressed  and  controlled  by  the 
computer  via  the  Multiprogrammer  backplane  or  can  be 
controlled  by  external  triggers  through  rear  card  edge 
connectors.  [Ref.  10:pp.  1-1,  2-1,  3-1,  Ref.  1:  p.  39] 

The  system  used  in  this  experiment  was  configured  with  one 
HP69759A  500  KHz  A/D  Converter  card,  one  HP69791A  High  Speed 
Memory  card  and  four  HP69792A  Memory  Expansion  cards. 

1.  500  KHZ  A/D  card 

The  500  KHz  A/D  card,  when  triggered,  converts  input 
voltages  into  12 -bit,  twos-complement  binary  values  within 
2.0  ns.  It  may  be  triggered  internally  by  program  command  or 
externally  by  a  lockout  input  at  the  edge  connector.  Digital 
output  is  available  at  the  edge  connector  at  a  rate  of  500 
KHz.  [Ref.  11] 

The  system  was  configured  to  receive  analog  data  from 
the  diode  array  and  send  digital  data  to  memory  at  the  edge 
connector.  Additionally,  it  was  configured  for  external 
triggering  from  the  A/D  trigger  circuit. 
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2.  High  Speed  Memory  and  Memory  Expansion  Cards 

The  High  Speed  Memory  card  can  store  up  to  65,536  16- 
bit  words  at  high  speed.  Data  can  be  written  to,  or  read 
from,  the  card  by  the  computer  via  the  backplane  or  an 
external  device  via  the  card  edge  connector.  Memory  Expansion 
cards  are  added  to  the  Memory  card  by  a  special  chaining 
cable.  Each  expansion  card  can  store  up  to  196,608  16-bit 
words.  With  four  expansion  cards,  the  total  memory  available 
was  851,968  16-bit  words.  [Ref.  12] 

Memory  was  programmed  to  operate  in  a  First  In/First 
Out  (FIFO)  mode.  This  allowed  the  cards  to  act  like  a  data 
buffer  with  data  written  to  sequential  memory  locations. 
[Ref.  12] 

Data  from  the  A/D  card  was  written  to  memory  via  the 
edge  connector.  Data  was  then  read  from  memory  to  the 
computer  hard  disk  by  the  computer  via  the  Multiprogrammer 
backplane . 
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C.  A/D  TRIGGER  CIRCUIT 
1 .  Purpose 

The  most  critical  part  of  the  Data  Acquisition  System 
was  the  A/D  trigger  circuit.  The  circuit  design  had  to  insure 
that  the  Diode  Array  Clock  signal  (500.0  KHz),  the  Diode  Array 
Blanking  Pulse  signal  (484.5  Hz),  and  the  "Take  Data"  command 
(200.0  Hz)  from  the  chopper  were  properly  combined.  The 
resultant  trigger  signal  was  required  to  be  a  group  of  1024 
pulses  (at  500  KHz)  which  were  synchronized  with  the  start  of 
a  diode  array  sweep  and  the  leading  edge  of  a  chopper  disc 
aperture.  Figure  4  is  a  diagram  showing  how  the  signals  were 
combined . 
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A/D  Trigger  Circuit 
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2.  Modifications 

During  previous  applications  a  Timer/Pacer  card  was 
programmed  as  an  input  to  the  A/D  trigger  signal,  and  a 
circuit  card  for  the  trigger  was  constructed  [Ref.  l:p.  41, 
Ref.  9:p.  42].  Since  the  Timer/Pacer  portion  was  deleted  and 
the  chopper  "Take  Data"  signal  (an  analog  voltage  signal  from 
a  photodiode)  was  added,  a  new  circuit  card  had  to  be  designed 
and  constructed,  as  previously  discussed. 

The  new  circuit  card  was  installed  and  successfully 
tested.  Figure  5  is  a  timing  diagram  of  the  A/D  trigger 
circuit.  A  detailed  drawing  of  the  new  circuit  card  is 
provided  in  Appendix  B. 

clock  1I!1||||||!||||||||||J||I!I|||||||||||U 

Blank  r  -y  y—y  -  y  y—y—  y—y - y— j 

blank  i _ i _ n._j_  n _ i _ i _ i _ o _ n _ u 

SKI  i _ i _ i _ n _ J 

trigger  _  mini  -  imu _ mum  miiitm _ 

Figure  5 

A/D  Trigger  Timing  Diagram 


21 


V.  SYSTEM  TEST  PROCEDURE  AMD  RESULTS 


A.  TEST  PROCEDURE 

1.  Initial  Alignment 

The  system  was  initially  constructed  with  the  light 
source  and  spectrograph  at  a  10.0  ft.  separation.  The  optics 
were  then  aligned  and  diode  array  voltages  checked  for  both 
the  reference  and  scene  beams.  The  variable  beam  splitter  was 
adjusted  to  provide  optimum  light  intensities  for  both  beams. 

It  became  readily  apparent  that  the  200  watt  light 
source  provided  insufficient  radiant  power  to  overcome  both 
being  passed  through  two  beam  splitters  and  the  losses  in  the 
fiber  optic  cable.  With  the  best  beam  splitter  adjustment  and 
system  optical  alignment  diode  array  peak  voltages  were  barely 
discernable  above  base  noise  levels.  The  data  acquisition 
system  was  unable  to  record  usable  readings  for  transmittance 
calculations . 

2.  Test  Configuration  and  Procedure 

Due  to  time  and  cost  constraints,  it  was  not  possible 
to  obtain  a  more  powerful  light  source.  The  system  was 
reconfigured  so  that  sufficient  diode  array  voltages  were 
obtained  to  prove  design  concept.  The  variable  beam  splitter 
was  removed  and  light  was  directed  through  the  chopper  into 
the  fiber  optic  cable.  The  beam  splitter  at  the  entrance  to 
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the  spectrograph  was  removed  and  the  output  of  the  fiber  optic 
cable  was  directed  into  the  spectrograph.  Reference  beam 
intensities  were  obtained  and  averaged.  It  was  then  assumed 
that  reference  beam  intensities  were  constant  for  the 
remainder  of  the  system  test. 

The  system  was  then  configured  to  obtain  test  data 
from  a  particle  source.  The  fiber  optic  cable  was  removed  and 
the  scene  beam  was  directed  through  the  chopper  and  directed 
onto  the  spectrograph  inlet  slit.  Additionally,  the 
spectrograph  was  moved  to  distance  of  approximately  4.0  ft. 
from  the  light  source  in  order  to  insure  sufficient  scene  beam 
intensity.  Calibration  (no  plume)  intensities  were  recorded 
prior  to  the  introduction  of  soot  particles.  The  smoke  cloud 
from  a  fuel-rich  oxy-acetylene  torch  flame  was  directed  into 
the  scene  beam  and  intensities  were  recorded. 

B.  TEST  RESULTS 

1 .  Data  Reduction 

The  first  step  in  data  reduction  was  to  obtain 
calculated  extinction  coefficient  (Q)  ratios  using  the  Mie- 
scattering  computer  code.  Q  was  determined  as  a  function  of 
d32,  o,  m,  and  wavelength.  Five  wavelengths  were  observed 
during  the  test  (X3  =  366.3,  X2  =  404.7,  X3  =  435.8, 

=  546.1,  and  X5  =  577.0  nm)  which  provided  ten  independent 
ln-transmittance  ratios. 
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Three  values  of  m  were  used  in  the  calculations.  The 
first,  m  =  1.8  -  0.3i,  is  similar  to  the  value  derived  by 
Kunitomo  and  Sato  for  various  soots  and  the  value  measured  by 
Chippett  and  Gray  for  propane  and  acetylene  soots .  The 

second,  m  *  1.95  -  0.66i,  is  the  complex  index  of  refraction 

of  carbon  measured  by  Senftleben  and  Benedict.  The  third, 
m  *  1.8  -  0.6i,  is  an  average  of  the  index  for  carbon  and  the 
index  for  propane  and  acetylene  soots  measured  by  Dalzell  and 
Sarofim.  Two  values  of  o  were  used,  1.5  and  2.0,  which  are 
typical  of  soot  particle  size  distributions.  [Ref.  4] 

The  measured  intensities  were  then  converted  into  ln- 
transmittance  ratios.  The  average  of  six  sweeps  of  the  diode 
array  were  used  to  determine  the  reference  beam  intensities 
for  five  observed  wavelengths.  The  intensities  for  the  scene 
beam,  both  with  and  without  particles,  were  similarly 

calculated  using  the  averages  from  two  sweeps  of  the  diode 
array.  The  transmittances  for  each  wavelength  were  then 

cai  ciliated  as  discussed  in  Chapter  II.  The  natural  logarithms 
were  computed  and  the  ten  independent  ln-transmittance  ratios 
(for  five  wavelengths)  were  obtained. 

The  next  step  in  data  reduction  was  to  look  for  rough 
trends  in  data  to  reduce  the  amount  of  required  regression 
analysis.  Values  for  d32  were  calculated  using  the  measured 
ln-transmittance  ratios  and  various  specified  values  of  m  and 
o  as  shown  in  Table  1.  The  refractive  index 
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m  *  1.8  -  0.6i  with  o  =  1.5  or  2.0  provided  a  poor  correlation 
since  six  of  the  ten  measured  ln-transmittance  ratios  were 
lower  than  the  minimum  calculated  extinction  coefficient 
ratios.  Additionally,  m  =  1.95  -  0.66i  with  o  =  2.0  yielded 
poor  results.  These  three  data  sets  were  eliminated  from 
further  consideration.  Calculations  of  mean  d}2  indicated  a 
likely  value  ranging  from  approximately  0.10  to  0.17  microns. 
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TABLE  1 


Calculated  d„  for  Specified  m  and  o  Values 


1  m  =  1.8  -  0.3 

S  m  =  1.95  -0.66 

!  m  x  1  8  -  0.6  1 

In  -  ratios 

19BE9 

ffSSSil 

S5SH1 

0.4884 

0.14 

0.06 

0.06 

a  a 

a  a 

a  a 

0.5865 

0.15 

0.09 

0.1 

a  a 

a  a 

a  a 

0  6106 

0  18 

mum 

0.125 

0.085 

0  12 

0  065 

0  6296 

0.14 

0.07 

0.09 

a  a 

a  a 

a  a 

0.7332 

0.21 

0.17 

0.155 

0  12 

0.15 

0  115 

0.7758 

0.135 

0.08 

0.09 

a  a 

a  a 

a  a 

0.7872 

0.22 

0.2 

0.155 

0  12 

0.16 

0  115 

0.7999 

a  a 

a  a 

a  a 

a  a 

a  a 

a  a 

0.8328 

0.09 

0.06 

0.06 

a  a 

#  a 

a  a 

0  9135 

0  2 

0.12 

0.12 

0.08 

0  12 

0  065 

msan  dlam. 

0.163 

0.108 

0.106 

0.101 

0.138 

0.00 

std  dev. 

0.04 

0.047 

0.036 

0  019 

0.018  0  025  | 

Note:  all  values  of  d32  are  in  microns 

**  indicates  ln-transmittance  ratio  below 
minimum  extinction  coefficient  ratio 


26 


Now  that  a  likely  range  of  values  for  d32  was  obtained, 
the  ratios  of  calculated  mean  extinction  coefficients  for 
assumed  values  of  m,  o,  and  d32  were  plotted  against  observed 
values  of  the  ln-transmittance  ratios.  If  the  ratios  had 
corresponded  exactly  with  each  other,  the  data  points  would 
have  formed  a  45  degree  line  passing  through  the  origin. 
Short  of  this,  the  best  fit  of  data  would  have  the  least 
scatter  about  this  line,  and  would  form  a  line  with  a  slope 
approximately  equal  to  one.  [Ref.  1] 

The  best  fit  of  data  was  found  by  applying  a  linear 
least-squares  fit  to  the  data.  The  regression  analysis 
function  in  Microsoft  Excel  was  ideal  for  this  purpose.  The 
correlation  coefficient,  r,  is  a  measure  of  the  relationship 
between  the  ln-transmittance  ratios  and  the  extinction 
coefficient  ratios.  The  square  of  the  correlation  coefficient 
expressed  as  a  percentage  was  used  to  measure  how  closely  the 
data  points  matched  the  regression  line.  The  r2  values 
closest  to  one  indicated  the  best  linear  correlation  between 
data.  The  x  -  coefficient  (slope)  of  the  regression  analysis 
indicated  how  closely  the  slope  of  the  least-squares  fit 
matched  a  45  degree  line.  [Refs.  1,  13] 

Figures  6,  7  and  8  are  plots  of  the  regression  lines 
for  m  =  1.8  -  0.3i  (o  =  1.5),  m  =  1.8  -  0.3i  (o  =  2.0),  and 
m  =  1.95  -  0.66i  (o  =  1.5),  respectively.  Particle  diameter 
(d32)  was  varied  from  0.09y,  to  0.1 8|a  in  each  case. 
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•  bar  ratios 
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2.  Data  Analysis 

As  previously  noted.  Table  1  data  indicated  that 
m  =  1.8  -  0.6i  with  a  =  1.5  or  2.0  and  m  =  1.95  -  0.66i  with 
o  =  2.0  provided  poor  correlation.  No  correlations  were 
obtainable  using  a  ratio  of  In  T(Xs=577.0n)  over 
In  T(X4=546 . Ip) .  This  apparently  resulted  because  the 
wavelengths  were  so  closely  spaced  [Ref.  5]. 

Table  2  lists  the  r2  and  x-coef ficient  values  obtained 
in  each  case.  The  best  r2  value  was  0.8839  for  d32  =  0.12p, 
m  =  1.8  -  0.3i,  and  a  =  2.0.  All  other  r2  values  were  within 
2.0  percent.  The  best  x  -  coefficient  value  was  0.995  for 
d32  =  0.15p,  m  =  1.8  -  0.3i,  and  o  =  1.5.  This  was 
significantly  closer  to  1.0  than  any  of  the  other  results, 
indicating  that  d32  is  probably  closer  to  0.15p  than  0.12p. 
The  poor  correlations  indicate  the  possible  existence  of  a  bi- 
modal  particle  distribution  in  which  the  second  mode  is  at  a 
larger  diameter  as  a  result  of  soot  agglomeration  in  the  fuel- 
rich  post-flame  region,  where  the  measurements  were  made. 

It  must  be  noted  that  the  system  recorded  many  more 
sweeps  of  the  diode  array  than  were  actually  used.  Better 
results  will  be  obtained  by  using  all  available  data  to 
calculate  measured  transmittances . 
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TABLE  2 


Summary  of  rJ  and  x  -  coefficient  Values 


m  s  1.8  -  0.3  sigma  *  1  5 

18-  0.3  sigma  =  2.0 

m*1.95'0  66  sigma  =  15  | 

diameter 

BUSHES 

x  -  coeHicient 

mm 

x  -  coeHicient 

ma 

x  -  coe  five  vent 

0  09 

0.8809 

1  0538 

0.8837 

0  9406 

0  8806 

0  9715 

0.12 

0  8759 

1.0486 

0.8839 

0  8789 

0  8753 

0  8999 

0.15 

0  8739 

0.995 

0.8838 

0  8057 

0  8705 

0  7783 

0.18 

0  8699 

0  8919 

0.8834 

0.7317 

0  8643 

0  6408 
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3 .  System  Performance 

With  the  exception  of  an  inadequate  light  source, 
system  performance  appeared  to  be  excellent.  The  A/D  trigger 
circuit  and  light  chopper  functioned  as  designed.  One 
shortcoming  in  system  software  was  noted.  If  the  data 
acquisition  program  was  actuated  during  an  A/D  trigger  pulse 
group,  then  the  recorded  data  was  obtained  out  of  sequence. 
This  problem  was  mitigated  by  slowing  the  chopper  rotation 
frequency,  and  thus  the  frequency  of  the  "take  data"  signal. 
This  resulted  in  longer  nulls  between  A/D  trigger  pulse  groups 
and,  therefore,  less  chance  of  actuating  the  data  acquisition 
program  during  a  trigger  group.  The  slower  chopper  rotation 
rate,  and  resulting  slower  data  acquisition  rate,  could  result 
in  degraded  performance  in  a  high  vibration  environment. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  system  design  concept  was  proven  valid.  The  goals  of 
obtaining  rapid  measurements  of  reference  and  scene  beam 
intensities  and  of  extending  the  distance  between  light  source 
and  spectrograph  were  not  achieved.  Installation  of  a 
stronger  light  source  will  make  these  goals  possible.  Poor 
correlation  of  the  data  obtained  from  a  fuel-rich 
oxy-acetylene  flame  was  probably  the  result  of  the  existence 
of  a  multimodal  particle  size  distribution. 

B.  RECOMMENDATIONS 

Recommendations  for  further  system  development  and 
applications  are: 

(1)  Obtain  a  more  powerful  light  source.  Minimum  power 
should  be  1000  watts. 

(2)  When  adequate  intensities  are  obtained  in  both  the 
scene  and  reference  beams,  test  the  system  on  a  particle 
source  with  known  characteristics. 
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(3)  Conduct  tests  in  environments  with  varying  amplitudes 
and  frequencies  of  vibration  to  determine  system  limitations. 

(4)  If  the  system  is  sensitive  to  vibrations  with 
amplitudes  and  frequencies  expected  in  the  field  environment, 
it  may  be  necessary  to  reduce  the  number  of  diodes  swept 
(since  only  approximately  10  of  the  1024  are  needed).  This 
will  result  in  a  increased  sweep  rate  (since 

TiWtp  =  #  diodes  in  array  /  500  Khz).  Additionally,  a  chopper 
capable  of  rotating  at  higher  speeds  will  provide  a  higher 
data  acquisition  rate,  since  a  higher  "take  data"  signal 
frequency  would  result.  The  data  acquisition  program 
and/or  the  A/D  trigger  circuit  would  require  modifications  to 
eliminate  the  problem  of  triggering  the  A/D  card  in  the  middle 
of  a  trigger  pulse  group. 

(5)  Obtain  the  fiber  optic  cable  and  lenses  required  to 
extend  transmission  distance  to  that  required  for  testing  over 
longer  distances  (~40  ft). 
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APPENDIX  A 

LIGHT  CHOPPER  DISC  DESIGN 


APPENDIX  B 
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DATA  ACQUISITION  PROGRAM 

!«*»»»»'  ACQUIRES  HULTIPLE  SHAMS  OF  DIODE  ARRAY  iiiimimimmmi! 
""'ii’sn  AND  STORES  JIITM  W  DIRK  mxiuimittii* 

BERT  I  WIST  II,  KFN  GPANAM, KELLY  HARRIS  l?0<  immimmm 
iHiimmHimiiimotmi'tHulitiiHnmmiiiHMiiiiiiiiuiiiiyi 

immim  REVISED  FOR  USE  WITH  5I0KIIZ  A/D  AMD  STORAGE  iimnnnnmi 
immim  OH  HARD  DRIW  BY  JOHN  VAHT-IU  197?  iimmimiiiH 
Inmnmii  REVISED  FOR  DUAL  BEAU  APPLICATION  smnmmxmimmii 
inmimum  BY  KEVIN  T  A  YEAR  1973  m  xiHxtiixiiiMHmxxinixi 
PHI  CU:BI?2),DIl:in?2)lZ(l :8I?2) 

LOCAL  lUCKOUl  7  1  DISAHU  S  f  RUN!  PANELS  OF  DAI  A  ACQUISITION  SYSTEM 

RFNOIE  70? 

WJIPUE  70?j‘AR'  !  ANALOG  RESFT  OPENS  All.  CMAMMFIS 

GOSUB  Initiate 
PRINT  USING  "R" 

PJIJJT  ^OATA  WILL  PE  SIORED  OH  HARD  DRIW!  WITH  FOtlOVIMO  fill  HANES: ■ 

TRim  •  CALIBRATION  — -  A  * 

PRINT  ■  PARTICI  F  SANP1.E  -----  B  * 

OUTPUT  70?:*Di:iI,2"  !  OPEN  SHUT  If  R 

PRINT  USING  7/" 

PRINT  ’TAKING  DATA* 

! WAIT  3 

G05UB  H"ltipreq  Hales  calibration  data 
G05UB  Perqe 
GOSUB  Sttredata 
PRINT  USING  '////////* 

PRINT  "  INTRODUCE  PARTICLES  IIIEN  PRESS  KEY  I  9  TO  TAH  DATA' 

ON  KEY  7  LABEL  ’  IAKE  DATA"  GOTO  1011 
Standby: GOTO  Standby 
ParM 

PR  INI  IJSIHG  "///’ 

PRINT  MAKING  PARTICI  F.  DATA" 

WAIT  .1 


GOSUB  Nnltiprnq 

GIISIW  Storpdata 
UUIPltl  70? ;  "AR ' 

LOIAI  7 
GOTO  3000 
Noltiprsq: 

Cdelay-0 
0UIFUT  707;  "MU  1,2" 

Cdnlay-Cdplayl  j 
If  Cdela»<5  THEN 
GOSUB  Initiate 
GOTO  1117 
END  IF 

IF  Par  1-0  THEN 
PRINT  USING  *///’ 

PRINT  'CALIPRAflOH  DATA  IS  STORED  ON  IDHORY  CARDS 
UAH  3 
END  IF 

IF  Part-#  THEN  GOTO  1570 
PRINT  USING  V//* 

PR  INI  'PARTICLE  DATA  IS  STURFD  OH  Iff  WIRY  CAPDS  • 
RETURN 


!  TAKE  PARTICI E  DATA 
TARTir.II  DATA  FILE  NAME 

•A/D  is  in  slot!  1 
•H?n  Card  2,3 

! liner /nacer  3 

•  Hen  F.ip  Cards  <,5,6,7 
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I A90  Stiredali: ! 

I  Ml  PRINT  USING  ■///’ 

1682  PRINT  ’MIA  IS  PEINC  1RAHSFERF0  FRON  THE  UVDS  TO  IMF  CONFUTED’ 

I6IJ  PRINT  USING  V 

1M5  PRINT  'WHEN  THE  TRANSFER  IS  COMPLETE  HIE  FIRST  5  FNIREES  OF  EACH  FILE* 
1686  PRINT  ’  WILL  PE  DISPLAYED’ 

IM7  OU1PUI  7?3;'WF2.3,76T’  IFIFU-IN  WITH  LOCKOUT 

1608  OUTPUT  723j’NI’;2;8l?2;’T’ 

160?  ENTER  7?38f.iZ<«> 

1612  IF  Par! )=1  TIEN  GOTO  1700 

1613  PRINT  'NEGATIVE  VOLTAGE  LOOP’ 

1615  FOR  1-1  TO  81 ?2 

1616  IF  ZUK32880  TIEN 

1617  C(I)~Z(1>*.IBI5 

1618  GOTO  1631 
161?  END  IF 

1621  c(ib<655J5-zu»*.oi05 
1639  NEXT  I 
1651  GFTTO  1782 
1  /OH  PirbParHl 

1711  PRINT  "NEGATIVE  VOI  PAGE  LOW 

1718  ir  Pir»( i  Turn  goto  i6i? 

1721  FOP  hi  TU  8172 
1738  IF  Z(  I  >(32888  THEM 
1710  D(lbZ(I)».RR8j 


1711  GOTO  1751 

1712  END  IF 

1758  0(I)”(655.15-Z(I>)*.I8I3 
1751  NEXT  I 
1768  rart=| 

1775  PRINT  USING  ’///’ 

1782  L«ckyhCIWK65> 

1703  ASSIGN  PTUh2  10  li-ikyl 
1701  OUTPUT  |Pa|h2!r.(U 
1785  IF  ParbR  fHER  GOIO  1770 
1787  Iicky8-CI1R8(66) 

1708  ASSIGN  IPUh2  TO  LockyO 
178?  OUTPUT  |Path2:0(»> 

1770  FUR  h65  io  65«Part 
17?l  IvckvhCHRt(l) 

1792  ASSIGN  fP»ih2  TO  l»cky1 
1 7? 3  PRINT  *PDAT  ’jCHRKI) 
I7?1  FOR  J=l  TO  5 

T7?5  ENTER  pPa»h?;N 

1776  TRINT  H 


1777  NEXT  J 


1798  NEXT  I 
17??  WAIT  3 


1881  OUTPUT  7?3;’Wr2.3,56T" 
1082  RETURN 


IFirO-m  WITHOUT  LOCKOUT 


1041 

1081 

1881 

1882 

1881 

1884 

1885 

1886 
1887 
1 

1H88 
ward  > 
1889 
1098 

1891 

1892 

1893 

1894 

1895 

1896 

1897 

1898 

1899 
1988 
1981 
1902 
1983 

RB 

MB 


3,256.2.3,56.I*!initialiie  the  ftede  register 
1  reference  2(leu  order  war tf ) 

1,31,2.2,8,1'ISet  reference  ?(hiuh  erder  word) 

1  ,?6|2.2,l/7777.T,!Set  reference  1  ( law  erder  word) 
1,27,2.2,17. l*!set  reference  llliinh  erder  uerd) 
1,22,2.2,1,  rlfteset  Read  FointerUeu  erder  uerd) 


l,23,2.2,l,T“!Reset  Read  Pain  ter  (hiijh  erder  uurd) 
l,?l,2.2,8,T*!Reset  the  liff ereere  counter  (lew  or 


order  uerd 


Initiate;! 

1HI1PIIT  71*3}  *UF,  2 
Mil  PUT  723; -Uf  ,2 
OlllPUf  7?3;*Mf'  2 
OUIPII1  723;  *UE  ,2 
(WlPllf  7?3;*Uf,2 
OtllPUl  723; ‘Uf  ,2 
OUlPlir  7?l;*Uf ,2 
(NIIPUI  723;*UF,2 

OUTPUT  7;'3;*WF,?.l,21,2.2,0,T*!Rese1  the  Difference  ceenterdiigh  erder 

milPtll  7?3;*UF, 2. 1,24,2.2,8.! 
ompill  721;"UF  2.1.25.2.2  8  1 
081PII1  723{*CC  2,lf*!Wa  t 
I’RKlI  USUtC  *3/* 

OUTPUT  723;‘UF ,2.3,561* 

OUTPUT  7?1;  “CY4:.?:  "7  * 

OUTPUI  723;‘UE1. 1,11, 1.2,81* 


OUTPUT 

OUTPUT 

0II1PUI 

OUTPUT 

0U1PIJ1 

OUTPUT 

ouipw 
no fpu  i 

OUTPUT 

OUTPUT 

KETURM 


7:»3;  -UF1 
7?3;*UF1 
7?3;*UF1 
723; “UF1 
723; "WFI 
723;*Ufl 
7?3;*UF1 
723;  *UT 
723;*0F1 
723;’UF1 , 


.1,12,1.2,81* 
.1  13,1.2  81* 
3,81* 

1,2, 1.2, IT* 
.1,1,1.2,01* 
.1,3,1.2,81’ 
.1,4, 1.2  IT* 
l, 5, 1.2,1. I" 
1,7.1. 2.lt* 
1,2<,|J,1T* 


IReset  the  Write  pointerUeu  order  uord) 
"IReset  the  Wfrite  pointer  (high  order  uerd) 
leer  card  instruction 

'SET  HFHOttY  CARD  TO  FIFO- 111  M)DE 
lACJlVAfE  CARD 

IDlSABtE  MUX  (Heed  to  enable  nui  if  ever 
Isuitrh  hack  to  I417A  connecter  cable) 
•DRW.  OUTPUT  FALSE 
!EIWI  EFN  FAT  SE 
! RESET  A/D  MODE  MTS 
!SFT  EX1ERMAI  TR1CCFR 
IDISAME  INTERNAL  IlfCCtR 
ITRICCFR  POT  ARTTY  SFlECT  I1IPUT  D1SAWE0 
TEXT  L0C1C  SENSE  AT.TIDE  HIGH 
1EXIERI1AL  LOCKOUT-INPUT  EMABlEO 
!WX7 IPLE  EXTERNAL  TRIGGER  RESPONSE 
•PROGRAM  CAIN  XlOlFOR  LSB  OF  .0005) 


1916  Purge:  ! 

!?.?0  FOR  1=65  TO  66 
193U  Leclyl-CHRFd) 

1948  P1RCE  Lucky! 

1951  HFXT  I 

I960  FUR  1=65  TO  66 

1978  I  »ckyi=CI1R!U) 

19U8  LREA1E  PDA?  I  uiky!,0192,  It 

1981  ASSIGN  8Patl.2  Til  Lucky! 

1982  NEXT  I 

1983  PRINT  IlfilllG  *///* 

1985  l*R  HIT  "HAD  TO  PIJffCF  AND  CRT  ATE  BOAT  FILES* 
2888  RETURN 
3081  END 


40 


LI8T  OP  REFERENCES 


1.  Vaughn,  John  K. ,  Measurement  of  Submicron  A L203  Particles 
in  Rocket  Plumes,  Master's  Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  December  1992,  pp  1-34. 

2.  Few,  J.,  Lewis,  J.,  Hornkohl,  J.,  Optical  Measurements  of 
Turbine  Engine  Exhaust  Particulates,  University  of 
Tennessee  Space  Institute-Center  for  Laser  Application, 
Report  No.  NAPC-PE-221C,  July  1991,  pp  1-8. 

3.  Few,  J.,  Lewis,  J.,  Hornkohl,  j..  Optical  Particle  Sizing 
in  the  Exhaust  of  an  A GT  1500  Combustor,  University  of 
Tennessee  Space  Institute-Center  for  Laser  Application, 
Report  No.  AC  DIV-PE-235C,  July  1992,  pp  1-7. 

4.  Cashdollar,  K.,  Lee,  C.,  Singer,  J.,  Three-Mavelength 
Light  Transmission  Technique  to  Measure  Smoke  Particle 
Size  and  Distribution,  Applied  Optics,  Vol.  18,  No.  11, 
June  1979,  pp  1763-1769. 

5.  Kim,  H. ,  Multiple-Wavelength  Transmission  Measurements  in 
Rocket  Motor  Plumes,  Master's  Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  September  1991,  pp  3-8,40. 

6.  Few,  J.,  Hornkohl,  J.,  Measurement  and  Predication  of  Jet 
Engine  Particulate  Effluents,  University  of  Tennessee 
Space  Institute-Center  for  Laser  Application,  Report  No. 
NAPC-PE-195C,  February  1990,  pp  1-11,  17-20. 

7.  Oriel  Corporation,  Light  Sources,  Monochronometers , 
Detection  Systems,  Catalog,  Vol.  II,  Statford, 
Connecticut,  1990,  pp.  82-94. 

8.  EG  &  G  Reticon,  Reticon  Linear  Diode  Arrays,  Application 
Notes,  No.  121,  Sunnyvale,  California. 

9.  Harris,  R.,  An  Apparatus  for  Sizing  Particulate  Matter  in 
Solid  Rocket  Motors,  Master's  Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  June  1984,  pp  40-41. 

10.  Hewlett-Packard  Corporation,  Data  Capture  with  the 
6942A/6944A  Multiprograimers ,  Application  Note  316-5, 
December  1984. 


41 


11.  Hewlett-Packard  Corporation,  H-P  500  Kilohertz  Analog-to- 
Digital  Converter  Card  Model  69759A,  Operating  Manual, 
June  1984,  p  1-2. 

12.  Hewlett-Packard  Corporation,  H-P  High  Speed  Memory  Card 
Model  69791A  and  Memory  Expansion  Card  Model  69792, 
Operating  Manual,  March  1984,  pp  1-1, 1-2. 

13.  Microsoft  Corporation,  Microsoft  Excel  User's  Guide  2, 
Redmond,  Washington,  1992,  pp  41-45. 


42 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  California  93943-5000 

3.  Department  Chairman,  Code  AA  1 

Department  of  Aeronautics  and  Astronautics 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

4.  Professor  David  W.  Netzer,  Code  AA/Nt  2 

Department  of  Aeronautics  and  Astronautics 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

5.  Professor  Richard  M.  Howard,  Code  AA/Ho  1 

Department  of  Aeronautics  and  Astronautics 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

6.  LCDR  Kevin  B.  Taylor,  USN  2 

438  E.  Fair  Harbor  Circle 

Houston,  Texas  77079 


43 


